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A distribution function for localized carriers, f(E, T) 
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is proposed by solving 



e (E-E a) /k B T +Ttr/ . 



a rate equation, in which, electrical carriers' generation, thermal escape, recapture and radiative 
recombination are taken into account. Based on this distribution function, a model is developed 
for luminescence from localized state ensemble with a Gaussian-type density of states. The model 
reproduces quantitatively all the anomalous temperature behaviors of localized state luminescence. 
It reduces to the well-known band-tail and luminescence quenching models under certain approxi- 
mations. 
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Carrier localization is a common phenomenon in many material systems such as semi- 
conductor alloys, quantum wells (QWs) and self-assembled quantum dots (QDs). It has 
profound effects on electrical and optical properties of the materials. It has been known 
for long that a number of anomalous temperature-dependent luminescence behaviors are 
related to carrier localization, including 1) the "S-shaped" temperature dependence of the 
luminescence peak positio n 1 ' 2 ^ 4 ! 5 and 2) a reduction of luminescence linewidth with increas- 
ing temperature. 5 ' 6 Till now, no theory based on microscopic viewpoint is present, which 
is capable of providing satisfactory explanations for all the anomalous behaviors of the 
luminescence. In this work, we propose a distribution function for localized carriers and 
develop a model for the luminescence of localized state ensemble (LSE). The model suggests 
that thermal redistribution of localized carriers within the localized states is the cause of the 
anomalies in the temperature dependence of the luminescence peak. It is further shown that 
the current model reduces to the well-known band-tail model 2 - at high temperature and to 
the luminescence quenching model of a two-level system^ when the distribution of localized 
states approaches a 5-function. The model is applied to quantitatively interpret the lumi- 
nescence data of two kinds of material systems, namely InGaAs/GaAs self-assembled QDs 
and InGaN/GaN QDs. In the former system, LSE naturally form due to size-distributed 
QDs while in the later, LSE can form due to the In-rich clusters with broad size distribu- 
tion. The physical implication of E a — E and its effect on the temperature dependence of 
luminescence peak position are revealed. 

For a system with localized electronic states having the density of states (DOS) p(E), 
the rate of change of carriers' population density N(E, T) in the state at energy E and 
temperature T is given by^ 



diV <f r > = G{E) + 2^ f AE) - "&Vj>-*>/** - (1) 

(it J\ Ti r T r 

where E a is the energetic position of a delocalized state to which the localized carriers ther- 
mally escape. The first term on the right, G(E), represents the rate of carrier generation, 
which is proportional to p(E) according to G(E) = k ■ p{E), 8 where k is a constant pro- 
portional to absorption cross section and the number of incident excitation photons. Note 
that the carrier generation includes both excitation of carriers instantly in localized states 
and capture of carriers those be excited outside. The second term on the right represents 



2 



the number of carriers re-captured by the localized states per unit time, in which 7 C is the 
recapture coefficient N' is the total number of carriers that are thermally activated away 
from the localized states as given by 

N , = r + °° N(E',T) e{E ,_ Ea)/kBTdE/ ^ (2) 

J-oo T tr 

in which l/r tr is the escape rate of the localized carrier. A = J™ p(E')dE' is the total 
number of localized states. The third term on the right hand side of Eq. gives the 
thermal escape rate of the localized carriers and the last term describes the de-population 
rate of the carriers due to radiative recombination. The term l/r r represents the rate of 
radiative recombination. r tr and r r are assumed to be the same for all the localized states. 
In Eq. the tunneling transfer of carriers^ between QDs is not taken into account because 
the non-resonant tunneling rate is negligibly small.^ More detailed discussion will be given 
later. The solution of Eq. ((Tj) under steady-state condition (dN/dt—0) is 



N(E, T) = A - n(E, T), (3a) 



where 



a _ KT tr /ou\ 

~ (1 - 7c) + (VwhA/A' 1 ' 

'^^ ^Ji^^ . (3c) 

In Eq. (|3b"jl. 6 = n{E f , T)dE' . Further expressing n(E, T) as n{E, T) = f(E, T) ■ p(E), 
then a distribution function can be derived as 

/( E ' T ) = e (E-E a )/k B T + Ttr / Tr - W 

Note that n(E,T) essentially describes the shape of the luminescence spectrum given by 
N(E,T)/r r , because A is a function of T only. Such a quantity has previously been 
used to describe the lineshape of phonon-assisted exciton luminescence peaks in some polar 
semiconductors . 11112 

Considering a general case of a LSE system with Gaussian-type DOS, which may result 
from, e.g., fluctuations of quantum dot size or alloy composition: 

p(E) = p e-^- E ^\ (5) 
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in above, po and Eq are the amplitude and peak energy position, while a is the standard 
deviation of the distribution. In Fig. [H the quantities of p(E)/p and (r tr /r r )/(_E, T) are 
plotted with respect to E for both E a > E and E a < E . As will be shown in the following, 
both cases exist in real physical systems as reported in the literature. 

As mentioned earlier, n(E, T) represents the "shape" of the luminescence spectrum and 
the luminescence peak position can be found by setting dn(E,T)/dE = 0. We found that 
when 

E = E -x- k B T, (6) 

n(E,T) reaches its maximum. The temperature-dependent coefficient x(T) can be obtained 
by numerically solving the following equation: 



xe x 



a 



x 



| e (E -E a )/k B T ^ ^ 

T tr 



It is not difficult to see that Eq. (JJJ has one and only one solution for < x < (a/k B T) 2 . 

It should be noted that the temperature dependence described by Eq. (jOJ) is only due to 
carriers' thermal redistribution within the localized states. It is known that the bandgap of 
an idealized semiconductor material without localized electronic states is itself temperature 
dependent, which is usually described by the Varshni empirical formula.— After taking into 
account this factor, the variation of the peak position of luminescence from LSE is then 
given by 

E(T) = E - ^ f - x ■ k B T, (8) 

where a is the Varshni parameter and G the Debye temperature of the material. For the 
cases where the thermal redistribution of localized carriers is dominant in the temperature 
range studied, the temperature dependence of the luminescence peak can be well described 
with Eq. 

At high temperature region an approximated solution of Eq. ([2]) is found to be (cr/A^T) 2 .— 
Eq. (jSJ) then becomes 

E(T) =Eq- — , (9) 

y J + T k B T' K ' 

which is simply the band-tail model proposed by Eliseev et al. 2 So, the widely adopted 

band-tail model can be viewed as an approximation of the current model under high T 

region. 



Full width at half maximum height (FWHM) is another parameter for a luminescence 
spectrum, which is also embedded in n(E, T). The FWHM T C (T) of n(E, T) can be obtained 
by numerically solving n(E, T) = n(E p k, T)/2. As will be shown below, the variation of r c (T) 
with temperature exhibits a "valley", i.e., T C (T) decreases first and then increases with 
raising temperature. Besides the variation in line width due to the thermal redistribution 
of carriers within LSE, the broadening due to phonon and impurity/imperfection scattering 
should be taken into consideration. The effective FWHM of luminescence peak is determined 
by making convolution of n(E,T) and a Lorentzian function, [4E 2 + (r + r p ^) 2 ] -1 . Here 
T is due to impurity/imperfection scattering and T p h = o XT + r )Lo/[e hu3L °^ kB ' T — 1] due to 
phonon scattering.— In the RRS model 15 , a a and 7^0 are the acoustic- and optical-phonon 
coupling strength, respectively. 

The integrated intensity of the luminescence spectrum is proportional to the total number 
of localized carriers, i.e., 

/+00 r+00 
N(E',T)dE' = A n(E',T)dE'. (10) 
00 J —00 

Utilizing an integral approximation^ 



■y _|_ e a(x+b) y _|_ g2.416sin6» 

where 9 =arctan(a/2.41), Eq. (|Tnj) can be derived as^ 



dx ~ n ... . (11) 



/(T)oc a + (l- 7c )-exp 



(E Q - E a ) + k B T ■ ln(r r /r tr 



(12) 



^(k B Tf + 2(a/2.41) 2 

For cr=0 (i.e., a ^-functional DOS), the above expression is reduced to the well-known model 
describing thermal quenching of luminescence for a two-level system.^ Indeed, for the case 
of a = 0, no thermal re-distribution takes place and the system becomes essentially an 
equivalent two-level system. This fact thus further validates the current model, which is 
more general. Note that for the two-level system, E a — Eq is in effect the thermal activation 
energy of the carriers. 

Figure El shows the calculated profiles for N(E, T) at different temperatures. The 
parameters used in the calculation are £0=1.185 eV, -E a =-Eo+0.073 eV, cr=13 meV, 
T~tr/ T r = 0- 027/250. From the figure, it is seen that the typical anomalies in the temperature- 
dependent luminescence are reproduced. Fig. HHa) plots the peak positions of the spectra as 



a function of temperature together with that predicted by Varshni empirical formula using 
the Varshni parameter a=0.48 meV/K and the Debye temperature 0=270 K. The sum of 
the two contributions is given by the solid curve, which is seen to agree excellently with the 
experimental data for an Ino.35Gao.65 As QD sample.— The QD density of the sample em- 
ployed in the present work is ~5x 10 10 cm -2 . The average distance between QDs is estimated 
as 20 nm. Other details of the sample have been previously described elsewhere.— The non- 
resonant tunneling rate of carriers between QDs, which is estimated to be ~10 5 s _1 using 
the Wentzel-Kramers-Brillouin approximation,— is not taken into account in the model. It 
is known that the tunneling rate depends weakly on temperature. However, the thermal 
escape rate increases exponentially with temperature. For example, at 50 K, the thermal 
escape rate of carriers occupying QDs with high energy levels (i.e., 58 meV below E a ) 
reaches to 10 7 s _1 . Therefore, the tunneling transfer is disregarded in the model developed 
for interpretation of temperature-dependence of luminescence of LSE. Figure H^b) presents 
the dependence of the FWHM on temperature, from which, it is seen that the reduction of 
FWHM in the luminescence spectra is mainly due to the effect of redistribution of localized 
carriers. The effect of phonon scattering is to broaden the spectra monotonously as the tem- 
perature increases, whereas carrier thermal re-distribution results in a dependence showing 
a valley as already mentioned earlier. The combination of the two effect leads to the anoma- 
lous dependence of FWHM on temperature shown by the solid curve in Fig. EKb). In the 
calculations, the values of r =5.0 meV, oa=10 /zeV/K, 71,0=18.8 meV and ^^0=36 meV 
were adopted . 1 ^ 2( | ) Finally, Figure Efc) gives the integrated intensity of the luminescence. It 
can be found that when the recapture coefficient 7 C =0.9 is taken, the calculated intensity 
agrees well with the experimental data. 

The results presented above is for E a — Eq > 0. The value of E a — Eq measures the mean 
thermal activation energy for the localized carriers. For the system of InGaAs self-assembled 
QDs investigated in this work, the delocalized state is seen to locate at an energy 73 meV 
above the central position of the localized states. The origin of such a delocalized state may 
lie on the presence of a wetting layer due to the S-K growth mode of InGaAs on GaAs?*^ 

For the case of E a — E Q < 0, calculated luminescence peak positions as a function of 
temperature is given in Fig. H] for a few selected values of E a — E . The other parameters 
used are also listed in the figure, which remain unchanged for the whole calculations. It is 
noted that the "S-shaped" temperature dependence of luminescence peaksi can be modeled 
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well. The solid squares stand for the experimental data for an InGaN/GaN QWs.— In 
fact, the experimental data for InGaN/GaN QWs samples reported by other groups 2,3 can 
be quantitatively interpreted within the whole temperature range using the current model. 
Further, for the partially ordered GalnP epilayer system grown on GaAs,— the observed 
anomalous temperature dependence of the anti-Stokes photoluminescence peak can also be 
reproduced with the current model.— 

Finally, we briefly discuss the physical meaning of E a . Like the Fermi level in the Fermi- 
Dirac distribution function, E a in the distribution function derived in the work gives a 
"marking" level below which all the localized states are occupied by carriers. Its relative 
position to E essentially determines anomalous temperature dependence of luminescent 
peak position. For the cases of InGaN/GaN QWs (E a < E ), E a may be the quasi- Fermi 
level of samples,— which depends on concentration of carriers optically/electrically injected 
and magnitude of the built-in electric field within the samples. It is obvious that the value 
of E a — E depends individually on sample. As shown in Fig. HI the magnitude of E a — E 
determines the details of peak position variation with temperature. Another point we want 
to mention is that if E a is far below Eq, i.e., Eo — E a > 8a, the model will be no longer valid. 
Under such condition, the luminescent density, jV ^ a ' T - > ~ 10~ 5 (photon number per second), 
could be too weak to observe when a typical value of r r =l nanosecond was considered. 

In conclusion, a model is developed, which quantitatively describes the temperature- 
dependence of the luminescence spectra from localized carriers. It reproduces almost all the 
anomalies observed for the luminescence of LSE. It is demonstrated that the two well known 
band-tail and luminescence quenching models are simply the approximations of the current 
model under certain limiting conditions. 
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Figure Captions 

FIG. 1. Normalized distribution function and Gaussian-type density of states for localized 
carriers. Note that there can be two cases: E a — Eq > and E a — Eq < 0. 

FIG. 2. Calculated population densities of localized carriers N(E, T), which resemble lumi- 
nescence spectra, as a function of energy and temperature for the case of E a — E > 0. 
The curves are shifted along vertical direction for clarity. 

FIG. 3. Calculated temperature dependence of luminescence peak position (a); FWHM 
(b); and integrated intensity (c) for the case of E a — E > using the parameters 
given for Fig. |21 The squares are experimental data and the solid lines are calculated 
using corresponding equations as denoted. 

FIG. 4. Calculated temperature dependence of luminescence peak position for the cases of 
E a — E < 0, depicting the "S-shape" dependence curve. The squares are experimental 
data, the solid lines are calculated by using Eqs. (JHJ) and (J7J) by setting different values 
of E a — E as denoted. 
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FIG. 4: of 4. Q. Li, et al. 
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